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ABSTRACT 

We report the first results of a program to measure accurate stellar velocity dispersions in the bulges of the host 
galaxies of active galactic nuclei (AGNs) for which accurate black hole (BH) masses have been determined via 
reverberation mapping. We find good agreement between BH masses obtained from reverberation mapping, and 
from the M, - a relation as defined by quiescent galaxies, indicating a common relationship between active and 
quiescent black holes and their larger-scale environments. 

Subject headings: galaxies: active — galaxies: Seyfert 



1. INTRODUCTION 

Two techniques have been successfully applied over the last 
few years to derive masses of supermassive black holes (BHs) 
in samples of galaxies. Kinematics of stars or gas on scales 
of < 1 - 10 pc yield estimates of M, based on Newton's laws 
under the assumption that the observed motions reflect the grav- 
itational influence of the BH. This technique has been applied 
to galaxies near enough {d < 100 Mpc) that the radius of influ- 
ence of the BH, r/, = GM^/crl, is likely to have been resolved on 
the angular scales accessible to ground- and space-based tele- 
scopes, 1" > 6* > 0". 1. The second technique, "reverberation 
mapping" (RM), makes use of broad emission-line variability to 
estimate BH masses in active galactic nuclei (AGNs). The size 
of the broad emission-line region is estimated from time delays 
between variations in the continuum (assumed to originate near 
the BH) and the response of the broad emission lines (Bland- 
ford & McKee 1982; Peterson 1993; Netzer & Peterson 1997). 
By combining the time delays with the emission-line widths, 
the BH mass then follows from the virial theorem. Black hole 
masses based on RM have been published for about three dozen 
AGNs (Wandel, Peterson & Malkan 1999; Kaspi et al. 2000). 

Both mass-estimation techniques have strengths and weak- 
nesses. Masses derived from RM are based on gas that is much 
closer to the BH, often r < 0.01 pc, than the stars or gas used 
in kinematical studies. At these separations, the gravitational 
force to which the gas responds is due almost entirely to the 
BH and there is little danger of a "false detection." However 
the precise relation between M, and the measured quantities is 
uncertain at levels of ^ 50% due primarily to uncertainties in 
the geometry and kinematics of the BLR and other unknowns 
(e.g.,Krolik 2001). 

In the case of BH masses derived from spatially resolved 
kinematics, formal accuracies can be high provided r/, has ac- 
tually been resolved; otherwise, estimates of M, depend sen- 
sitively on assumptions about the character of the motion near 
the BH. Here the danger is that the inferred BH mass is actually 



a measure of the combined mass of the BH and the stars in the 
resolved region, r ^ ri,; modeling techniques applied to such 
data can substantially overestimate M, by ascribing too much 
of the gravitational force in the nucleus to the BH and too little 
to the stars. 

In fact, there have been claims of systematic differences be- 
tween RM and kinematical BH masses, in the sense that RM 
masses at a given bulge luminosity, Mb, are a factor ~ 20 
smaller than kinematical masses (Wandel 1999). The discrep- 
ancy has most often been ascribed to systematic errors in the 
RM masses (Richstone et al. 1998; Faber 1999; Ho 1999). Un- 
fortunately, it is extremely difficult to make a direct compar- 
ison between dynamical and reverberation-based BH masses. 
In AGNs, the variable nonstellar continuum and broad emis- 
sion lines that enable RM tend to swamp the stellar features at 
r < ri, which are needed to study the stellar dynamics. At the 
present time, there is no galaxy for which the BH mass has been 
measured using both techniques. 

A possible solution to this problem emerged with the discov- 
ery of the M, - a relation, a tight empirical correlation between 
M, and the stellar velocity dispersion a (Ferrarese & Merritt 
2000, hereafter FMOO; Gebhardt et al. 2000a, hereafter GOOa). 
The latter is defined on scales much larger than r/, and hence is 
easily accessible to ground-based instruments. Since stellar ve- 
locity dispersions can also be measured in active galaxy hosts, it 
is possible to determine whether the M, - a relation for AGNs 
is similar to that of quiescent galaxies. Preliminary attempts 
at such a comparison have been published (Gebhardt et al. 
2000b [hereafter GOOb]; Merritt & Ferrarese 2001a [hereafter 
MFOla]) based on published AGN bulge velocity dispersions 
(principally from Nelson & Whittle 1995, hereafter NW95). 
These studies reveal a general consistency between the AGN 
data and the M, - a relation for normal galaxies, confirming 
that the apparent discrepancies between the kinematic and RM 
masses were due mostly to systematic overestimates in the BH 
masses derived from stellar kinematic measurements (MFOla, 
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Merritt & Ferrarese 2001b). Unfortunately, the uncertainties in 
both the AGN BH masses and stellar velocity dispersions used 
by GOOb and MFOla are too large to provide a truly critical test. 
The reason for the large uncertainties in M, is quite simple: 
the largest source of random errors in AGN BH masses is the 
light-travel time delay for the emission lines, which determines 
the size of the line-emitting region. The fractional error in this 
quantity is largest for those sources that have the smallest line- 
emitting regions, since the uncertainties are determined in large 
part by the time resolution of the monitoring program. The 
smallest line-emitting regions are found in the least-luminous 
objects (i.e., R oc lP-^~^ ''-^ Kaspi et al. 2000). Generally speak- 
ing, lower-luminosity AGNs have not been monitored any more 
closely than higher-luminosity AGNs and thus the fractional er- 
rors in their BH masses are relatively large, in some cases as 
large as 100%. The low-luminosity AGNs are the very galaxies 
for which bulge velocity dispersions have been measured, and 
thus do not yet provide a conclusive argument for consistency 
between the M, - a relations in AGNs and normal galaxies. 

In order to better compare the M, - a relations for AGNs 
and quiescent galaxies, we are undertaking a program to mea- 
sure bulge velocity dispersions of the host galaxies of AGNs 
with formally well-determined reverberation-based BH masses. 
In this contribution, we describe the first results from this pro- 
gram. 

2. OBSERVATIONS AND DATA REDUCTION 

2. 1 . Sample Selection 

Black hole masses based on reverberation studies have been 
published for 34 AGNs (Wandel, Peterson & Malkan 1999; 
Kaspi et al. 2000). There are both systematic and random un- 
certainties in these measurements: systematic errors arise pri- 
marily from the largely unknown geometry and kinematics of 
the line-emitting region, and the random errors are due mostly 
to uncertainties in the size of the line-emitting region, as noted 
earlier, and the uncertainty in the translation of the observed 
FWHM to the 3-D velocity dispersion. The typical level of 
random uncertainty in the reverberation-based BH masses is 
^ 50%. We focus our attention on those AGNs with better- 
determined BH masses (i.e., smaller random errors). These 
provide the most meaningful comparison with normal galax- 
ies: the scatter in the M, - a relation for normal galaxies is 
^ 30%, while the intrinsic scatter in the relation appears to be 
negligible (FMOO). 

2.2. Observations 

We have obtained bulge velocity dispersions from the Call 
triplet lines at rest wavelengths 8498, 8542, and 8662 A. These 
lines are in the spectral region where the AGN contribution is 
minimized (e.g., NW95). 

The objects observed in this study are listed in Table 1. Col- 
umn (1) gives the common name of the object, and its redshift is 
given in column (2). Column (3) gives the reverberation-based 
virial mass of the central BH in each object, and the published 
source of the virial mass is given in column (4). 

The observations reported here were obtained on UTC 2001 
April 7-9 on the Mayall 4-m Telescope of the National Opti- 
cal Astronomy Observatories on Kitt Peak. The observations 
were made with the Ritchey-Chretien Spectrograph with the 
BL380 grating (1200 Hnes per mm, blazed at 9000 A, cen- 
tered at approximately 8910 A), and an RG610 blocking filter 
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The detector employed was the T2KB CCD, a back-illuminated 
2048 X 2048 device windowed to 440 x 2048 pixels. The slit 
width was set to 2". The resulting spectra have a resolution 
R K. 5000, and cover the range 8 150-9660 A. To allow accu- 
rate sky-subtraction, a long slit (~ 5') was used. Column (5) 
of Table 1 gives the spectrograph position angle at which the 
observations were made; when such information was available, 
we aligned the slit along the major axis of the galaxy. Column 
(6) gives the total integration time on each source. 

Each observation of a galaxy was bracketed with quartz-lamp 
spectra for flat-field correction, and He-Ne-Ar spectra for wave- 
length calibration. A spectrum of the twilight sky was obtained 
to determine a nightly cross-dispersion illumination correction. 
We also obtained spectra of 14 late-type giant stars (G8 III — 
K6 III) that were used as templates and radial-velocity stan- 
dards in the analysis described below. 

2.3. Data Reduction 

The raw long-slit spectra were reduced by flat-field, bias, 
and slit-illumination corrections following the procedures de- 
scribed by Pogge (1992). We used the XVlSTA^ package for 
all reductions. Curvature of the long-slit spectra in the dis- 
persion and cross-dispersion directions was removed by tracing 
stars and bright, isolated night-sky emission lines in the spectra 
and rectifying the long-slit spectra. In this spectral region, the 
T2KB CCD is subject to fringing at the 10% level longward of 
^9000 A, which we reduced by using the bracketing flat-field 
exposures. In general, residual fringing was less than 0.1% in 
the spectral regions of interest. We combined the longslit spec- 
tra in groups of three 20-min exposures to remove cosmic -ray 
events by sigma clipping with a simple CCD noise model (gain 
= 1.9 e", (Tread-out = 5e"). The nuclear spectra were extracted 
from each group within a 2" x 4" aperture centered on the nu- 
cleus. For the fainter objects, we used a 4".8 extraction aperture 
to improve the signal-to-noise ratio. Sky spectra were extracted 
from flanking regions along the slit and subtracted from the nu- 
clear spectra. We then combined all of the nuclear spectra into 
a final spectrum. This group-wise extraction of spectra taken 
within at most an hour of each other reduced the effects of vari- 
ability in the atmospheric OH airglow lines. The fully reduced 
spectra of all objects are shown in Fig. 1. 

We point out that the M, - a relation was defined by FMOO 
using central velocity dispersions corrected to an aperture of 
size rg/8, with the effective radius of the galaxy. GOOa, on 
the other hand, used spatially averaged, rms, line-of-sight ve- 
locity dispersion within r^. MFOla demonstrated that there is 
no systematic difference between the two definitions of a in this 
sample. The bulge effective radius is not known for any of our 
target galaxies, but the bulge is enclosed within the aperture 
used for the spectral extraction. On this basis, in §3 we treat 
our measured dispersions on the same basis as the aperture- 
corrected dispersions tabulated by MFOla. 

2.4. Data Analysis 

Two independent estimates of the stellar velocity dispersion 
were obtained from the Call triplet. The first employed the 
Fourier Correlation Quotient (FCQ) method described by Ben- 
der (1990) and Bender, Saglia, and Gerhard (1994). The sec- 
ond was based on the Maximum Penalized Likelihood (MPL) 
method developed by Merritt (1997). FCQ and MPL differ in 
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Table 1 
AGN Sample Observed 



Galaxy 
(1) 


Redshift 

z 

(2) 


M. 

(3) 


Source^ 
(4) 


Position 
Angle 
(5) 


Exposure Time 
(sec) 
(6) 


cr(FCQ) 

(kms-i) 
(7) 


(t(MPL) 
(kms-^) 
(8) 


NGC 4051 

NGC 4151 
NGC 5548 
Mrk79 
MrkllO 
Mrk817 


0.0024 
0.0033 
0.017 
0.022 
0.036 
0.031 


14+010 

^ -^"-0.70 

5.9 ±2.5 

10.2!^:^ 

0.77!«:i 

■-'•■-'^-0.86 


1 

2 
3 
2 
2 
2 


135° 
135° 
110° 
56° 
110° 
90° 


3600 
3600 
7200 
7200 
15600 
12000 


80±4 
93±5 
183 ±10 

130±9 
86±5 

142 ±6 


80±3 
85±5 
180 ±6 
120 ±8 
95±8 
140 ±8 



^Data source for masses: 1: Peterson et al. (2000); 2: Kaspi et al. (2000); 3: Peterson & Wandel 
(2000) 
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Fig. 1 . — Calibrated spectra for the six galaxies in our program. An arbitrary 
additive constant has been added to the fluxes on the vertical axis. The dotted 
line represents the best fit returned by the FCQ method. 



important ways (e.g., Joseph et al. 2001); using both allows 
us to assess the impact of systematics on the measured values. 
MPL was applied to the fully calibrated and continuum nor- 
malized spectra. FCQ incorporates a continuum fitting routine, 
thus was applied to the spectra before normalization. The ve- 
locity dispersion measurements resulting from FCQ and MPL 
are Usted in columns (7) and (8) of Table 1, respectively. The 
errors associated with a are the formal random uncertainties 
returned by the FCQ and MPL codes, but there are likely unac- 
counted systematic uncertainties. Repeating the analysis using 
different stellar templates, or changing the parameters for the 
continuum fits in FCQ show that a more realistic error estimate 
is of order 15%. 

3. RESULTS AND DISCUSSION 

No velocity dispersions have been previously published for 
Mrk79, MrkllO and Mrk817. NW95 hsted velocity disper- 
sions from Call and/or Mgb absorption Unes for six galaxies 
with reverberation-based BH masses, including NGC4151 and 
NGC 5548. Unfortunately, NW95 warn that only for two of 
these galaxies, NGC 4051 and 3C120, are the measured dis- 
persions deemed accurate. For NGC 4051, they quote a = 
88 ± 13 kms~\ in good agreement with our result. For 3C120, 
which has a RM mass of 3.0^} 4 x 1O''M0 (Wandel, Peterson 
& Malkan 1999; Kaspi et al. 2000) NW95 fist a = 162 ±20 
km s~' (from Smith, Heckman, & IlUngworth 1990). For lack of 
an alternative, the NW95 velocity dispersions for all six galax- 
ies were adopted by GOOb and MFOla, in spite of their large 
associated uncertainties. 

In Fig. 2, we show the relationship between BH mass and 
bulge velocity dispersion (FCQ values with 15% errors as- 
sumed) for the six galaxies listed in Table 1, along with the 
same relationship for quiescent or weakly active galaxies (full 
references for both M, and a for the latter can be found in 
MFOla). For completeness, we also plot, as small open circles, 
the sample of galaxies for which values of M, were tabulated by 
GOOa based on dynamical modeling of Hubble Space Telescope 
data, although we refrain from using these data in the following 
discussion since the modeling and analysis of the data that led 
to those M, measurements have yet to be published. 

While secure conclusions are premature given the small sam- 
ple size, several points are worth mentioning. As noted by 
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GOOb and MFOla based on lower-quality data than presented 
in this paper, the current evidence is that reverberation-based 
masses are not systematically underestimated as suggested by, 
e.g., Richstone et al. (1998), Faber (1999) and Ho (1999). We 
stress that the data presented in this paper are unique for having 
high-quality measurements of both M, and g. Our data support 
the conclusion that dynamical and reverberation-based masses 
are generally consistent, and that the two methods can poten- 
tially yield BH mass estimates of comparable precision. This is 
of particular relevance. Reverberation mapping is intrinsically 
unbiased with respect to BH mass, provided the galaxies can be 
monitored at closely spaced time intervals. While dynamical 
methods rely on the ability to spatially resolve the region dom- 
inated by the BH gravitational potential, RM samples a region 
which is per se unresolvable. It follows that an aggressive mon- 
itoring campaign of AGNs is the only viable method to probe 
the low-mass (M, less than a few million solar masses) end of 
the M, - a relation. It is also apparently as good a method as 
traditional dynamical studies to probe the higher-mass regime, 
although the obvious drawback is that it is perforce limited to 
galaxies with Type 1 AGN, which constitute only ~ 1% of the 
general galaxy population. Probing the low mass end of the M, 
- (J relation is of particular interest since the slope and scatter 
of the relation have important implications for both hierarchi- 
cal models of galaxy formation (Haehnelt, Natarajan, & Rees 
1998; Silk & Rees 1998; Haehnelt & Kauffmann 2000) and the 
effect of mergers on the subsequent evolution (Cattaneo et al. 
1999), including the coalescence of the binary BHs which are 
expected to form as a result of galaxy mergers (Milosavljevic 
& Merritt 2001). 

Furthermore, RM can probe galaxies at high redshift and 
with a wide range of nuclear activity, opening an avenue to 
explore possible dependences of the M, - a relation on red- 
shift and activity level. Regarding the latter point, we note that 
NGC 4051, the galaxy with the smallest M, in our sample, is 
a narrow-line Seyfert 1 (NLSl) galaxy. Mrk 110 also approxi- 
mately meets the traditional narrow-line Seyfert 1 criterion that 
FWHM(H/3) < 2000 km s-'. The BH masses of the NLSls are 
consistent with those of the other AGNs, and all are consistent 
with the masses inferred from the M, - a relation defined by the 
quiescent galaxies. This result is at odds with the recent find- 
ings of Mathur et al. (2001), who claim that NLSls lie system- 
atically below the M, - a relation defined by Seyfert 1 galaxies. 
One paradigm for NLSls is that, compared to AGNs with sim- 
ilar non-thermal luminosity, they are powered by lower mass 
BHs accreting at a relatively higher accretion rate. A compet- 
ing explanation is that they differ from other AGNs only in in- 
clination; the small widths of the broad lines and concomitantly 
low reverberation masses are due to viewing disk-like sources 
at low inclination (i.e., nearly face-on). The consistency of the 
BH mass and bulge velocity dispersion of NGC 405 1 , in partic- 
ular, with the M, - a relation for other galaxies strongly favors 
the low-mass/high accretion rate interpretation for the nature of 
NLSls. 

A linear fit to the sample of 18 galaxies — 12 galaxies from 
FMOO plus the six reverberation-mapped galaxies — using re- 
gression with bivariate errors and intrinsic scatter (Akritas & 
Bershady 1996; see also discussion in MFOla) gives a slope of 




(T (km ') 

Fig. 2. — Black hole mass versus central velocity dispersion a of the host 
elliptical galaxy or bulge. Small filled circles and errorbars represent quiescent 
or weakly active galaxies with dynamical measurements of M, . Open circles 
represent unpublished dynamical BH masses listed by GOOa. Large solid tri- 
angles are the reverberation-mapped galaxies presented in this paper, while 
open triangles are reverberation-mapped galaxies for which cr is available in 
the literature. The solid line is the best linear fit to the M, - a relation as pub- 
lished by MFOla for the kinematical masses only (filled circles), with slope 
4.81 it 0.55. The dotted line is the best fit to the six RM galaxies presented in 
this paper (fiUed triangles), and has a slope 5.3 ± 1.6. 



4.50 ±0.34, consistent with the slope of 4.81 ±0.55 derived 
by MFOla using only the kinematical mass estimates. To con- 
clude, we note that three of the galaxies in Fig. 2, NGC 4151, 
Mrk 79, and Mrk 1 10, seem to lie slightly above the "fiducial" 
M, - a relation defined by MFOla; i.e., these galaxies have 
black holes larger than expected from their bulge velocity dis- 
persions. This is unlikely to be attributable to systematic er- 
rors in a, since the quality of the fit returned by both FCQ and 
MPL for these galaxies, while not as good as for NGC 405 1 and 
NGC 5548, is comparable to that of Mrk 817, which Ues on the 
fiducial relation. There also seems to be no correlation between 
residuals and distance^. A more detailed discussion is prema- 
ture at this point: the small discrepancy could be due to small 
systematic errors in the RM masses, or, perhaps more interest- 
ing, may indicate that the current level of activity has allowed 
for significant growth of the BH mass. These results underscore 
the importance of further testing the M, - a relation by mea- 
suring a for a larger sample of reverberation-mapped galaxies 
over a large BH mass range. 

We are grateful to NOAO for observing time on the 4-m 
telescope, and especially to H. Halbedel and B. Gillespie for 
capable assistance at the telescope. LF acknowledges NASA 
LTSA grant NAG5-8693. This research has made use of the 
NASA/IR\C Extragalactic Database (NED), which is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



*Note that reverberation-mapping masses are independent of distance, therefore the scatter cannot be due to errors in the distance determination to these galaxies. 
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